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ABSTRACT 

In proto-neutron stars with strong magnetic fields, the neutrino-nucleon 
scattering/absorption cross sections depend on the direction of neutrino momentum 
with respect to the magnetic field axis, a manifestation of parity violation in weak 
interactions. We study the deleptonization and thermal cooling (via neutrino 
emission) of proto-neutron stars in the presence of such asymmetric neutrino opacities. 
Significant asymmetry in neutrino emission is obtained due to multiple neutrino-nucleon 
scatterings. For an ordered magnetic field threading the neutron star interior, the 
fractional asymmetry in neutrino emission is about 0.006(5/10^^ G), corresponding 
to a pulsar kick velocity of about 200 (B/W^'^ G) kms^^ for a total radiated neutrino 
energy of 3 x 10^^ erg. 



Subject headings: stars: neutron — pulsars: general — supernovae: general — dense 
matter — magnetic fields — radiation transfer 



1. Introduction 



Recent analyses of pulsar proper motion (Lyne & Lorimer 1994; Lorimer et al. 1997; Hansen 
& Phinney 1997; Cordes & Chernoff 1998), observation of pulsar bow shocks (Cordes et al. 1993), 
and studies of pulsar-supernova remnant associations (Frail et al. 1994) indicate that supernova 
explosions are asymmetric and the neutron stars receive large kick velocities at birth (250-500 
kms~^ on average, but possibly with a significant population having velocities greater than 
1000 kms^^). Compelling evidence for supernova asymmetry also comes from the detection of 
geodetic precession in PSR 1913-1-16 (Cordes et al. 1990) and orbital plane precession in the PSR 
J0045-7319/B star binary (Kaspi et al. 1996). In addition, evolutionary considerations of neutron 
star binary population also imply the existence of pulsar kicks (e.g., Fryer et al. 1998). 
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Two classes of mechanisms for the natal kicks of pulsars have been suggested. The first class 
relies on local hydrodynamic instabilities in the proto-neutron stars (e.g., Burrows et al. 1995; 
Burrows & Hayes 1996; Janka &: Miiller 1996) as well as global asymmetric perturbations seeded 
in the presupernova cores (Goldreich et al. 1996). The asymmetries in matter and temperature 
distributions lead to asymmetric explosion and/or asymmetric neutrino emission, although the 
magnitude of the resulting kick velocity is still unclear. 

In this letter, we focus on the second class of models in which the large pulsar velocities arise 
from asymmetric neutrino emission induced by strong magnetic fields. Although a number of 
papers have been written on this subject (Chugai 1984; Dorofeev et al. 1985; Binovatyi-Kogan 
1993; Vilenkin 1995; Horowitz & Piekarewicz 1997; Horowitz &: Li 1997), they are unsatisfactory for 
a number of reasons: they either failed to identify the most relevant neutrino emission/interaction 
processes or the relevant physical conditions in proto-neutron stars, or stopped at estimating the 
magnetic field effects on neutrino opacities, and none of them treated neutrino transport in any 
systematic manner. As a result, their estimates of the kick velocity that may be produced are 
unreliable. We have carried out systematic studies of neutrino transport in strong magnetic fields 
(Lai & Qian 1998a, b) in order to pin down the magnitude of the magnetic field effects on neutrino 
transport and pulsar kicks. Here we report our most important finding on the macroscopic 
consequence of asymmetric neutrino opacities due to parity violation in weak interactions. More 
details can be found in our forthcoming publication (Lai & Qian 1998a, hereafter LQ). 



2. Parity Violation and Asymmetric Neutrino Opacities 

In the presence of an external magnetic field, the differential opacity for neutrino-nucleon 
scattering, z^(e,/x,r) + -/V — > V(e,n,T) + ^ can be written in the form: 



die ' \ 1 

1 + einfi • B + eout^^' • B + const x (fi • fi') , (1) 



where and Q,' are the unit momentum vectors for the incoming and outgoing neutrinos, 

(sc) 

respectively, B is the unit vector along the magnetic field direction, and is the total scattering 
opacity without magnetic field (the subscript "E" implies that the opacity depends on the neutrino 
energy E). The presence of the B-depedent terms in Eq. (|l]) implies that there is a preferred 
direction in the opacity. This is a manifestation of parity violation in weak interactions. The 
coefficients ein and Cout are related to the nucleon polarization P by 

^ _^j^ ca{ca-cv) _ ^ ca{ca + cv) 

where cy = —1/2, ca = —1.15/2 for neutron and cy = 1/2 — 2sin^6'iy = 0.035, ca = 1.37/2 
for proton (Raffelt & Seckel 1995). For nondegenerate nucleons, the spin polarization P is 
P = fJ-mB/kT = 3.15 X 10^^(7m-Bi4 (10 MeV/T), where fXm is the nucleon magnetic moment 
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{9m = —1-913 for neutron and 2.793 for proton), and i?i4 is the field strength in units of 10^^ G. 
In the degenerate regime, we have P = 3fimB/i2nN) = 4.73 x IQ-^gmBu (10 MeV//Li7v), where 
/^AT is the nucleon chemical potential (or Fermi energy). A general expression for P is given in LQ. 
In proto-neutron stars, one should average Eq. (|^) over composition (neutrons and protons). It is 
also important to note that the differential opacity for i'e,fi,T scattering on nucleons is obtained 
by interchanging ein and eout in Eq- (0) as a result of the crossing symmetry of the leading order 
matrix elements. 

The neutrino absorption opacity is also asymmetric, i.e., k^^^^^ = k^^^^\i + e^fJ • B), with 
Ca ~ P- However, in the bulk interior of the neutron star, this asymmetry is exactly canceled by 
the asymmetry in emission (see §4). One also expects asymmetry in v -\- e ^ v + e. However, 
for relativistic degenerate electrons, the v-e opacity is much smaller than the u-N opacity, while 
the electron polarization P^, — 2 x 10~^i?i4(200MeV//ie)^ is of the same order as P for nucleons. 
Therefore the contribution from v-e scattering to asymmetric neutrino transport is small compared 
with that from iz-N scattering. 



3. Toy Problem: Effect of Multiple Scatterings 

The asymmetry coefficients (ein and eout) in the opacity are extremely small even at i?i4 = 10. 
However, cumulative effect due to multiple scatterings can enhance the asymmetry in neutrino 
emission. To demonstrate this effect, we consider a simple toy problem in which neutrinos leak 
out of a slab (with infinite extent perpendicular to the z-axis which coincides with the magnetic 
axis) via biased diffusion, i.e., the neutrino scattering probability has the form cx (1 + ecos6'), 
with e <C 1 and 6 being the angle between the z-axis and the direction of neutrino propagation 
after each scattering. The evolution of the neutrino number density n{z, t) is governed by the 
Fokker-Planck equation: 

^ = icAV2n-lcV.(6nB), (3) 

where A is the (constant) mean free path of the neutrino. The surfaces of the slab are located 
at z = ±Ars (so 2ts is the total optical depth of the slab). Suppose we start with a uniform 
distribution n = no at t = (so that there is no initial asymmetry). As neutrinos gradually 
diffuse out, the profile n{z, t) becomes increasingly asymmetric, and the flux asymmetry AF/F 
grows (where F is the sum of the fluxes from both surfaces and AF is the difference). For 
t ^ Tg{X/c), the "cooling" front reaches the center of the slab, and AF/F saturates to constant 
value, AF/F ~ €Ts/2 (valid for er^ <C 1). This constant value can be derived analytically from a 
quasi-steady-state analysis. The fractional asymmetry in the total neutrino radiation energy is 
0.34 erg. 

In proto-neutron stars, the neutrino optical depth r is large. A rough estimate can be made 
by setting the diffusion time ~ 3i?^/ (Ac) {R is the neutron star radius) to 1 s (the spread in the 
arrival time of SN 1987A neutrinos was 10s), which gives r ~ 10^. This would imply a large 
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enhancement in neutrino emission asymmetry. Several complications arise when one applies the 
result of the toy problem to real neutron stars: As shown in §4, the asymmetry in neutrino flux is 
determined by the coefficient Cout — and the local neutrino energy distribution function. Due 
to the crossing symmetry, we have (eout — f^m)u = — (cout — £in)u- Because ^'^(t) and i^^(t) have 
the same local energy distribution function, the "drift fluxes" (see Eq. |^) associated with I'^^r) 
and j^^(r) exactly cancel. The "drift fluxes" of z^e and do not cancel because a newly-formed 
neutron star is lepton-rich, i.e., it contains more Ve than Pg; and the typical energy of i^e is 
also different from Pg. As the neutron star is deleptonized, the difference between and z^e 
becomes smaller, and we expect the asymmetry in the total neutrino flux to gradually diminish. 
Another issue which is not clear from the toy problem concerns the role of neutrino absorption: 
for electron-type neutrinos, absorption is more important than scattering (by a factor of a few) 
in the deleptonization phase. Does the neutrino absorption completely wipe out the asymmetry 
associated with multiple scatterings? 



4. Thermal Evolution of Proto-Neutron Stars with Asymmetric Neutrino Opacities 

We now describe our calculation of the thermal evolution of proto-neutron stars in the 
presence of the asymmetric neutrino opacities discussed in §2. We start from the general transport 
equation for the spectral intensity Ie = /£;(r, f2,t) of a given neutrino species: 




jr ^^-JE)a 

(4) 

where I'^ = lE{r,Q' ,t), /e = {2tt)^Ie/E'^ is the neutrino occupation number, and /^^^ is the 
Fermi-Dirac distribution function for neutrinos: 

^E U ) - ^^^E - (2^)3 e(E-^.)/T + 1 

{U^^^ is the corresponding energy density), with matter temperature T and neutrino chemical 
potential In Eq. (|^ the absorption opacity k^^^^^ already includes the effect of stimulated 
absorption of neutrinos. Multiplying Eq. (^ by and then integrating over ft, we obtain the first 
order moment equation: 



(abs) (sc) 

VUe + ea^^c {Ui^""^ -Ue)B + (eout - ^m)^ cUe{1- Ie) B, (6) 



o p ^ <j ^ E 



where Ue, F^; is the spectral energy density and flux, and = k"^^ + Hi~^^^ is the total opacity. 
The neutrino number flux is simply given by S^; = Fe/E. Equation (|6|) reveals that in addition to 
the usual "diffusive" flux (the flrst term), there is a "drift" flux (the B-dependent terms) along the 
magnetic fleld direction due to the asymmetric neutrino opacities. Note that when Ue = U^^^ 
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(i.e., when neutrinos are in thermal equilibrium with matter), the asymmetric absorption opacity 
does not contribute to the drift flux, since in this case asymmetric emission exactly cancels 
asymmetric absorption. In the following we shall focus on the bulk interior of the star (below the 
layer where neutrinos decouple from the matter) , where U e for i^e has the Fermi-Dirac form with 
chemical potential (ly^, while and Vr have zero chemical potential. 



At) 



For a given neutrino species i, the total opacity can be written in the form 



^(abs) 



_(abs) 



'^Ei ~ ^Ei ' ^ ^Ei"' ~ i^/^o) {^0 is a fiducial energy). In general, because of nucleon 
degeneracy and/or neutrino degeneracy, there can be additional i?-dependence in the prefactor 
Ki (Burrows & Lattimer 1986). In these cases, we evaluate Ki at an appropriate mean neutrino 
energy. If we make a further assumption that fe and Pe have the same form of opacities (Burrows 
& Lattimer 1986), then an analytical expression for the total energy flux (integrated over neutrino 
energy) associated with z^e and can be obtained: 



Fu. + Fr 



■V 




B, 



(7) 



where we have defined rj^,^ = fi^^/T, (i^ 
Eq. ^ we have used the relation = 
number flux is given by 



Ku^^ /ku}, and e 



and er 



= (cout - ein)t/,. In deriving 
Similarly, the net electron lepton 
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We shall group f^, u^, Ur and Vr together and denote them as f^. With opacity ke = Ku^iE / Eq)'^ , 
the total energy flux carried by these neutrinos is given by 



cEj 
18 p K^^ 



-VT^. 



(9) 



As noted earlier, the drift flux associated with p-type neutrinos is zero. The equations governing 
the thermal evolution of a proto-neutron star can be written as 



dt 
dU 

'di 



n 



— = -V-F 



V • (F,, 



+ + F, 



(10) 
(11) 



Here Yl = riL/n = + Yy = {Y^- — Y^+) + {Yy^ — YpJ is the lepton number fraction (n is 
the baryon number density), and U is the internal energy density of the medium (a sum of 
contributions from neutrinos, electrons, positrons, photons, and nucleons). Note that inside the 
star, /3-equilibrium holds to a good approximation, thus we have (1,^^ + pn = P-p + l^e — 1-293 MeV 
(lin, Pp do not include rest mass). For simplicity, we have neglected gravitational contraction and 
have adopted a Newtonian treatment. 
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We assume a uniform magnetic field in the z-direction throughout the neutron star. When 
the asymmetry is small, we can expand T{r,t) and YL{r,t) as 

T(r, t) = To(r, t) + AT(r, t) cos 6, YLir, t) = ¥^0(7; t) + AyL(r, t) cos 6, (12) 

and similarly for other relevant quantities. The angular dependence can then be factored out. We 
solve the evolutionary equations numerically using an explicit finite difference code. 

Our model neutron star has a mass M = I.SSMq and a radius R = 11km, with a density 
profile resembling a F = 3 polytrope (the central density is pc = 8 x 10^^ gcm~^). The initial 
temperature and Yj, profiles are shown in Fig. 1. In the stellar interior, we have = 0.28, 
Yj, = 0.366. The temperature is 8 MeV in the inner core (M ^ O.TM©) and reaches 22 MeV in the 
outer core due to shock heating. These initial conditions are typical of those found in core collapse 
simulations. The initial chemical potentials of various particles at the center are //^ = 97.7 MeV, 
fjLp = 52 MeV, He = 313 MeV, and jiy = 266 MeV. We assume that there is no asymmetry at 
t = 0. We adopt the zero boundary condition T = at the stellar surface. Since the neutrino 
flux asymmetry mainly accumulates in the interior of the star, our crude treatment of the surface 
boundary condition is adequate. 

Figure 1 depicts the Tq, Ylo profiles and the asymmetric perturbations (AT and AY^,) as 
functions of time for a magnetic field of strength S14 = 5. The evolution of the Tq and Y^q 
profiles is similar to that found by Burrow & Lattimcr (1986): The deleptonization phase (t ^ 10 
s) is dominated by the lepton number diffusion. In this early phase a heat wave moves inward, 
raising the core temperature. This is followed by a thermal cooling phase in which the core 
temperature decreases gradually due to energy diffusion. The asymmetries in the temperature and 
lepton number profiles increase during the deleptonization phase, reaching a maximum at t ~ 8 
s (This behavior is analogous to that found in the toy problem discussed in §3). This growth of 
asymmetry is also reflected by the increase in the ratio \/S.Lz/L\ (where L is the luminosity, and 
I^Lz/c is the net z-momentum carried away by the neutrinos per unit time) as shown in Fig. 2(b). 
As the neutron star becomes more lepton-depleted, the difference (in number and energy) between 
Ve and Ve diminishes, and the asymmetry gradually declines. Fig. 2(c) shows the asymmetry in 
the integrated neutrino luminosity [£{t) = J^dtL^t)] as a function of time. The asymmetry in the 
total radiated neutrino energy is a = \A£z/S\ ~ 0.006 5i4. This result is smaller than suggested 
by the toy problem (§3) mainly because the energy carried away by electron neutrinos during 
the deleptonization phase is only about 10% of the total neutrino energy released. In addition, 
geometric factors and cancellation of asymmetries associated with neutrons and protons tend to 
reduce the total asymmetry in neutrino emission. 

5. Conclusion 

In this paper, we have identified the most efficient way of generating natal pulsar kicks based 
on asymmetric neutrino emission induced by magnetic fields. The key point of this mechanism is 
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parity violation in weak interactions. The resulting pulsar kick velocity is 

(13) 

where -Etot is the total energy released by neutrinos (of all species) from the proto-neutron star, 
and (Bz) is the (averaged) ordered component of the magnetic field in the star. In reality, rotation 
(if it is misaligned with the magnetic axis) tends to reduce the kick velocity. We also note that the 
asymmetric 1/^ and emission from the proto-neutron star would give rise to asymmetric neutrino 
heating behind the stalled shock in the delayed supernova mechanism. This asymmetric heating 
could lead to asymmetric supernova explosions. 

Another related (but distinct) kick mechanism relies on the fact that Uf. and absorption 
opacities near the neutrinosphere depend on the local magnetic field strength (due to quantization 
of e~ and energy levels). If the magnitude of magnetic field in the north pole is different from 
that in the south pole, then asymmetric neutrino flux can be generated. We have shown in Lai & 
Qian (1998b) that this kick mechanism requires a much larger field strength than the mechanism 
considered in this letter. 

It is plausible that magnetic fields stronger than 10^^ G can be generated in proto-neutron 
stars (Thompson & Duncan 1993), and several pieces of evidence (albeit tentative) have been 
suggested in the literature (e.g., Vasisht & Gotthelf 1997). While hydrodynamically-driven 
kick mechanisms may still be viable, (depending on the magnitude of asymmetry seeded in the 
presupernova core), the possibility that a significant number of pulsars have kick velocities greater 
than 1000 km s~^ (e.g., Cordes & Chernoff 1998) may well require large magnetic fields (10^^ G) 
to be present in proto-neutron stars. 

We thank Sterl Phinney for useful discussion at the early stage of this work. D.L. also thanks 
Ira Wasserman for discussion. This work was started while D.L. held the Richard C. Tolman 
Fellowship at Caltech. Additional support was provided by NASA grant NAG 5-2756 and NSF 
grant AST-9417371. Y.Q. is supported by the David W. Morrisroe Fellowship at Caltech. 
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Fig. 1. — Thermal evolution of a proto-neutron star. The distributions of temperature and lepton 
number fraction (the left panels) and the asymmetric perturbations (the right panels) are shown 
as functions of time for a uniform magnetic field of strength Bu = 5. 
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Fig. 2. — The time evolution of (a) total neutrino luminosity L, (b) fractional asymmetry in 
luminosity AL^/L, and (c) fractional asymmetry in the integrated luminosity a, corresponding to 
the evolution depicted in Fig. 1. 



